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We systematically re-examine the extraction of the masses and couplings of the 1 hybrid, four-quark and molecule mesons 
from QCD spectral sum rules (QSSR). To NLO for the perturbative and power corrections, the hybrid mass is Mh ~ 1.81(6) 
GeV and Mh < 2.2(2) GeV from the positivity of the spectral function. In the same way , but to LO, the four-quark state 
mass is M4q = 1.70(4) GeV and M45 < 2.4(1) GeV, while the molecule mass is about 1.3(1) GeV. The observed 7ri(1400) and 
7ri(1600) might be explained by a two-component mixing with the set of input masses (1.2 ^ 1.3 ; 1.70 ~ 1.74) GeV and with 
a mixing angle 6 ~ — (11.7± 2.2)°, which slightly favours a molecule/four-quark mixing, and which eventually suggests that the 
7ri(2015) is mostly an hybrid meson. Isospin and non-exotic partners of the previous states and some of their radial excitations 
are also expected to be found in the energy region around 2 GeV. Further tests of this phenomenological scenario are required. 



1. Introduction 

Light hybrid mesons with the exotic quantum number 
1 ^ and with characterisitc decays into p7r, bin and 
rj'TT have been expected to be easily found in hadronic 
experiments. There are indeed two good experimental 
candidates 7ri(1400), 7ri(1600) T^, but their masses are 
lower than the lattice predictions in the quenched ap- 
proximation [2\ and with two dynamical quarks (2q) [3] : 

Mh\ quenched 

~ 1.9 GeV , Af^bq ^ 2.2(2) GeV . (1) 

First studied in QCD in |4l5j . lowest order (LO) results 
from QCD spectral sum rules (QSSR) [7] a la SVZ [5] 
have been obtained in [619110] : 



Mff|Lo-(1.6-2.1) GeV 



(2) 



which, after the inclusion of radiative corrections (partly 
obtained in ^llj and completed in [12 ) lead to the upper 
bound and estimate in units of GeV: 

A^^Inlo < (L9 - 2.0) , Mff Inlo = (1-6 ~ 1.7) . (3) 

Though the QSSR result tends to favour the hybrid as- 
signement of the 7ri(1600) meson, as quoted in [I], it is 
important to understand the discrepancy of the above 
QSSR results with the ones from the lattice and to a 
lesser extent with the one in [13] where it is claimed that 
the inclusion of the operator anomalous dimension ob- 
tained in [11] decreases substantially the previous result 
in Eq. ^ down to 1.2 GeV, though a such effect disap- 
pear to LO, and a similar one in the tensor meson has 
lead to negligible corrections [14]. One can also notice 
that [12] use as input the value of the four-quark con- 
densate from vacuum saturation which has been shown 
from different channels and different groups to be vio- 
lated by a factor 2-3 [7ll5ll6ll7ll8TTgj . In the following, 
we revisit the existing NLO results by paying attention 
on the influence of different QCD input parameters used 
in the analysis. We shall also check the effect of the op- 
erator anomalous dimension and study the new effect of 



the tachyonic gluon mass introduced in [11120121] but 

overlooked in [12]. For a more robust estimate, we shall 
consider results both from finite energy (FESR) and 
Laplace (LSR) sum rules. We complement our stud- 
ies by reexamining the predictions obtained for the 1 
light four-quark J2 and molecule _23J states and con- 
clude the paper by a phenomenological discussion on the 
possible nature of the observed 7ri(1400), 7ri(1600) and 
7ri(2015) mesons. 

2. QCD spectral sum rules (QSSR) 
Description of the method 

Since its discovery in 1979 ^8^, QSSR has proved to be a 
powerful method for understanding the hadronic proper- 
ties in terms of the fundamental QCD parameters such 
as the QCD coupling a^, the (running) quark masses 
and the quark and/or gluon QCD vacuum condensates 
[7] . In practice (like also the lattice) , one starts the anal- 
ysis from the two-point correlator (standard notations): 



(fx 



= -(<?^V-9^9'^)n« (g^) 



+g^^<;^(g2), 



(4) 



built from the hadronic local currents 0^^^{x): 



(5) 



which select the specific quantum numbers of the hy- 
brid mesons; A and V refer respectively to the vector 
and axial-vector currents. The invariant II'^^^ and 11*^°^ 
refer to the spin one and zero mesons. One exploits, in 
the sum rule approaches, the analyticity property of the 
correlator which obeys the well-known Kallen-Lehmann 
dispersion relation: 
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n: 



(1,0) 
V/A 



dt 



t — — ie 



- im n('^°) 



(6) 
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where ... represent subtraction terms which are poly- 
nomials in the q^-variable. In this way, the sum rule 
expresses in a clear way the duality between the inte- 
gral involving the spectral function lmllyj^{t) (which 
can be measured experimentally), and the full correla- 
tor U'^fjil'^)- The latt er can be calculated directly in 
QCD in the Euclidean space-time using perturbation 
theory (provided that —q'^ + m? (m being the running 
quark mass) is much greater than A^), and the Wil- 
son expansion in terms of the increasing dimensions of 
the quark and/or gluon condensates which simulate the 
non-perturbative effects of QCD. 

The SVZ expansion and beyond 

Using the Operator Product Expansion (OPE) [S], the 
two-point correlator reads for m = 0: 



n: 



V/A ' 



D=0,2, 



dimO—D 



where v is an arbitrary scale that separates the long- and 
short-distance dynamics; C are the Wilson coefficients 
calculable in perturbative QCD by means of Feynman 
diagrams techniques; {0{v)) are the quark and/or gluon 
condensates of dimension D. In the massless quark 
limit, one may expect the absence of terms of dimen- 
sion 2 due to gauge invariance. However, it has been 
emphasized recently [21] that the resummation of the 
large order terms of the perturbative series can be mim- 
iced by the effect of a tachyonic gluon mass A which 
generates an extra 1 /q^ term not present in the original 
OPE. This short distance mass has been estimated from 
the e+e^ data 18 20J and pion sum rule [50] to be: 

— ~ -(0.07 ± 0.03) GeV^ (7) 

TT 

In addition to Eq. ([7]) , the strengths of the vacuum con- 
densates having dimensions D < 6 are also under good 
control, namely: 

• (a,.G2) ~ (6.8 ± 1.3)10-2 GeV^ from sum rules of 
e^e~ ^ / = 1 hadrons [18j and heavy quarkonia 
[26124128] : 

. gii^Gi^) = g(^Aaa^''G^» - 2 x (0.8 ± 
0.1) GeV2(V'V'), from the baryons [29I19| and the 
heavy-light mesons [3D] systems; 

• as{i)il))^ ~ (4.5 ±0.3) X lO"'^ GeV^ from e+e" 
1=1 hadrons [TH] and r-decay data [T7], where 
a deviation from the vacuum saturation estimate 
has been noticed from different studies [7] and [151 

• g^{G^) ~ (1.2±0.2) GeY^iasG^) from dilute gaz 
instantons [31] and lattice calculations [32] . 

In the numerics, we shall use the value of the QCD scale: 

A3 = (353 ± 15) MeV , (8) 

deduced recently to 4-loops from the value of as{Mr) = 
0.3249(80) from r-decay [l7]. 



Spectral function 

In the absence of the complete data, the spectral func- 
tion is often parametrized using the "nai've" duality 
ansatz: 



1 



ImH 



(1,0) 
V/A 



Ml) 



"QCD continuum" x ^(t - tc) , (9) 



which has been tested ^ using e+e" and r-decay data, 
to give a good description of the spectral integral in the 
sum rule analysis even in the case of the broad a and 
states [33]; fn (an analogue to /^) is the hadron's 
coupling to the current; while tc is the QCD continuum's 
threshold. 

Sum rules and optimization procedure 

Among the different sum rules discussed in the literature 
[7], we shall be concerned with the following Laplace 
sum rule (LSR) and its ratio [8|25|26| : 

£^o)(r) = / dtrexp(-tr) -Imn|;;°)(t) , 



7^„(T) 



dr 



\ogCn , [n > 0) 



(10) 



The advantage of the Laplace sum rules with respect to 
the previous dispersion relation is the presence of the 
exponential weight factor which enhances the contribu- 
tion of the lowest resonance and low-energy region ac- 
cessible experimentally. For the QCD side, this pro- 
cedure has eliminated the ambiguity carried by sub- 
traction constants, arbitrary polynomial in q^, and has 
improved the convergence of the OPE by the presence 
of the factorial dumping factor for each condensates of 
given dimensions. The ratio of the sum rules is a useful 
quantity to work with, in the determination of the res- 
onance mass, as it is equal to the meson mass squared, 
in the usual duality ansatz parametrization. As one can 
notice, there are "a priori" two free external parame- 
ters (r, tc) in the analysis. The optimized result will 
be (in principle) insensitive to their variations. In some 
cases, the tc-stability is not reached due to the too naive 
parametrization of the spectral function. In order to re- 
store the tc-stability of the results one can fix the tc- 
values by the help of FESR (local duality) [27116] : 



Rn{tc) = 



1 

I" 

/J^citr+iilm<;°)(t) 



dt t" - imn[^/°V) 



/*^dtt«iimn|;;°)(t) 



M 



H 



(11) 



The results discussed in the next section will satisfy 
these stability criteria. 

3. The hybrid tv^ro-point function in QCD 

A QCD analysis of the two-point function have been 
done in the past by different groups [415] . where (unfor- 
tunately) the non-trivial QCD expressions were wrong 



3 



leading to some controversial predictions [?]• In this 
paper, we extend the analysis by taking into account 
the non-trivial correction and the effect of the new 
1/q^ term not taken into account into the SVZ expan- 
sion. The corrected QCD expressions of the correlator 
are given in I6l9l7j to lowest order of perturbative QCD 
but including the contributions of the condensates of 
dimensions lower than or equal to six. The new terms 
appearing in the OPE are presented in the following!^: 

• The perturbative QCD expression including the 
NLO radiative corrections reads ITT]: 



• To leading order in as, the contributions of the 
dimension-six gluon and mixed condensates read 



tt(0)| 



1 



487r2g2 
11 as 1 



00 _ 



18 ^ g2"^5(V^GV')log-^ , (16) 



where one can notice the miraculous cancellation 
of the log-coefficient of the dimension-six conden- 



sates for n 



(1) 



imn 



(1) I 



TTnn(°) I 

iiiiiiy^^ \pert 



607r2 



TT 



121 
16 



—n (— ^\ 1 — 



1207r2 
35 



1 + ^ 

TT 



1997 167 



432 360 



--^ log^ 
36 6 / ^ t 



(12) 



• The anomalous dimension of the current can be 
easily deduced to be [TT] : 



7 



71 



16 



(13) 



• The lowest order correction due to the (short dis- 
tance) tachyonic gluon mass reads: 



• The four-quark condensate contributions includ- 
ing radiative corrections read for n flavours |12j : 



lly 1 6 



1 16 

— —asiHA'^ 
9 



1 /91 
18 V"6" 



5n, 



Ols_ 
TT 



1/11 \a 



(17) 



• The contribution of the four-quark condensate in 
the (pseudo)scalar channels vanishes to leading or- 
der in as and starts at order a'^{ip'ip)'^ ■ In the 
scalar channel, it reads [12]: 



(0), 
V l6 



1 16 



asiH-'Y 



14 

y 



2 



as_ 

TT 



1/53 \as. (?2 

12 12+"^ y^°^-^ 



(18) 



Imn 



Imn 



(1) 

V/A 

(0) 

V/A 



{t)x = 



a, 35 , 



607r2 4 



as 15 



1207r2 2 



(14) 



• The contributions of the dimension-four conden- 
sates reads in the limit = to LO |6|9|7j and 
to NLO [12,: 



Imn 



V |4 



1 

9^ 

8q 

9~TT 



t 



log- 



imn^^^U 



1 

6tt 

SOs 
9 TT 



as{G^ 
t 



log- 



1 _ 145 as 

72 TT 



8asm{'ipTp) 



209 as 



8asm{ipip) 



(15) 



where as = as/TT and (tp^) = {ipuipu) — {fpdtpd) 



^ We neglect some possible mixings with operators containing 
more 7-matrices like g'tpi^a'y^cr^xxlijG'^^ , which is expected to be 
small. 



4. Upper bound on A//f(l +) from LSR 

Using the positivity of the spectral function, we can de- 
duce an upper bound on AIh from 7?.„(r) : ri = 0, 1. 
We show the result in Fig. [T]for the 1 ^ channel. The 
result is given in the domain limited by the two full (red) 
curves from TZq (t) and by the two green curves (dashed) 
from TZi{t) spanned by the two extremal values of the 
set of QCD parameters given in Table [1] 



Table 1 

Value of the upper bound on Mh in GeV at the minima of 'R.q{t) and 
at its intersection with 7^i(t) for two extremal values of the Sets of 
Power Corrections in units of GcV*^ {d is their dimension). 



Power Corrections 


(a,/7r)A2 




a^iiptpf X 10'' 




Set 1 


-0.04 


0.055 


4.8 


< 2.0 


Set 2 


-0.10 


0.081 


4.2 


< 2.4 



One can note that the prediction increases by 100 MeV 
each when |A2| and as(G2) decrease from their central 



1] 










/ 
/ 
















Figure 1. a) Domain of upper bound on M/f in GeV versus tiie LSR 
variable r in GeV~^ for two extremal values of the Sets of Power Cor- 
rections in Table [Tl from ■R-o(t) (red: continuous) and ?^i(r) (green: 
dashed); b) Effect of radiative corrections on the central value of the 
bound of Mh from 'R.o(r) red: : LO (dashed-dotted); NLO con- 
densate+LO perturbative (dashed) ; Total NLO (continuous). Green 
curves are from ?^i(r). 



with 7i the anomalous dimension in Eq. (jl3p and — /3i = 
1/2(11 — 2n//3) the first coefBcicnt of the /3 function for 
Uf flavours. 




value to the one allowed in the range given in Section 
2, while it increases by 25 MeV when 10^q;s(V'V')^ in- 
creases from 4.5 to 4.8. The effects of the variation of A 
and Mg within the range in Section 2 are invisible. At 
the stability points (in the sum rule variable r) of the 
continuous (red) curves which are also the intersections 
with the dashed (green) curves, one can deduce: 



Mh(1"+) < (2.2 ±0.2) GeV 



(19) 



This value confirms and improves previous results in 
|9I11I12| . The result (2.2 ± 0.2) GeV from lattice with 
two dynamical quarks though still consistent with this 
bound is on its boarder. We show in Fig [T] the effects 
of radiative corrections on the results. One can see that 
the ones due to the condensates give larger effects and 
increase the bound from 2.0 (LO) to 2.25 GeV, while the 
ones due to the perturbative terms decrease the mass 
prediction by 0.05 MeV, where the a^-correction and 
the one induced by the anomalous dimension act in the 
opposite directions. The ratio T^i is amost unaffected 
by the radiative corrections [dotted (green) curves on 
top of each others]. 

5. QSSR predictions of and Mh for 1"+ 

Following Ref. [13], we introduce the RGI coupling 
defined as: 



(logi//A) 



7i/-/3i 



(20) 



Figure 2. a) Domain of FESR predictions versus the QCD contin- 
uum threshold using the set of power correction in Tabled for jn 
from Alo (red: continuous) and A^i (black: dashed (on top of the red 
curve)) and A^2 (green: dotted); b) The same as for //j but for Mh'- 
from Ro{tc) (red: continuous); from Ri(tc) (green: dotted). 

Using FESR, we show the predictions for the decay con- 
stant in Fig. [2^) and the mass of the 1 in in Fig. [2)d). 
Using the stability criterion on the variation of tc, and 
the intersection of the predictions of different moments, 
and the Sets of Power Corrections in Table [1] one ob- 
tains the optimal value for the 1 



Jh = 9.6(1.4) MeV : = (4.4 ' 
Mh = 1.80(10) GeV : = (2.8 



5.2) GeV^ 
' 5.0) GcV^. (21) 



We show in Fig. [3] the predictions of the LSR in the 
1 ^ channel using the central values of the QCD input 
parameters and using two values of tc = 4.6 (beginning 
of the T-stability) and 5 GeV^ fixed from the previous 
FESR results. The prediction for the decay constant 
is not conclusive as there is not a stability in the LSR 
variable r. However, on can see from Fig. [3] that the 
different predictions interact for : 



Jh = (10 - 22) MeV , 



(22) 



which is consistent with the previous FESR result but 
less accurate. Therefore, we consider as a final predic- 
tion the one from FESR. For the mass prediction, the 
LSR gives: 

Mh - 1-81(3)(7) GeV : = (4.6 - 5.0) GeV^ (23) 

where the 1st and 2nd errors come respectively from tc 
and the QCD input parameters. We take the (naive 



Oti VtZ 0*3 0.4 a>e Dt4 Q.e o.e i 




Fi gUre 3. Predictions of the LSR in the 1 hybrid channel using 
the central values of the QCD input parameters and two values of tc — 
4.6 (beginning of the r-stability) and 5 GeV^ fixed from the previous 
FESR results: a) RG invariant coupling constant fn from Co (red: 
continuous), from £i (black: dashed-dotted (on top of the red curve)), 
from Ci (green: dashed); b) M/f from 7^o(t) (red: continuous) and 
from 7^1 (t) (green: dotted). 



Figure 4. Predictions of the four-quark state mass in GeV using 
the current in Eq. (1261 1 : a) LSR: domain spanned by the upper bound 
versus r using the largest range spanned by the QCD parameters given 
after Eq. (□ ( red: continuous curve); Mass using the central values 
of QCD parameters for tc — 3 GeV^ (green: dashed), — 4 GeV^ 
(black: dotted-dashed), ^ 5 GeV^ (blue: dashed), b) FESR versus 
tc'. region spanned by the largest range of QCD parameters. 



arithmetic) average of the FESR and LSR results and 
we take the quadratic average of the errors. Then, we 
obtain: 



(Mh) = 1.81(6) GeV 



(24) 



6. Comparison with existing theoretical results 



The central value of the result in Eq. ([M]) has increased 
by 100 - 200 MeV compared to the previous QSSR re- 
sults quoted in Eq. ([3]). The difference with the one in 
Ref. [9111] is due to the inclusion of radiative correc- 
tions for both perturbative and power corrections here, 
which increase the mass prediction by about 200 MeV. 
The difference with Ref. [Hj is the non-inclusion of the 
tachyonic gluon mass and the use of factorization for 
the four-quark condensate in [12]. In [13], a value of 1.2 
GeV for the mass has been obtained H, which the author 
attributes to be due to the hybrid operator anomalous 
dimension. In this paper, we show that the effect of per- 
turbative radiative corrections including the one due to 
the anomalous dimension only decreases the mass pre- 
dictions by 50 MeV. This result is (a priori) expected 
where it is easy to show that perturbative radiative cor- 
rections tend to cancel in the ratio of moments used 
to extract the meson mass. A similar explicit exam- 



^Also, in |13| an operator having high number of 7 matrices has 
been considered in the Ught-cone gauge but its correspondence in a 
covariant gauge is not quite transparent. The current looks to a be 
a tensor current while from the definition of the matrix element, 
the 1 ^ hybrid meson contributes through its longitudinal (i.e. 
spin zero) part. 



pie has been studied for the case of 2++gg current [14]. 
The present result is slightly lower than the lattice re- 
sults |2l3j quoted in Eq. ([1]). As the as corrections 
are reasonnably small, we do not see in the QCD side 
any potential contributions which can restore the dis- 
crepancy between QSSR and lattice results. From the 
phenomenological side, a more involved parametrization 
of the spectral function might help, but in many known 
channels, the usual duality ansatz : one resonance 4- 
QCD continuum describes accurately the spectral func- 
tion at the tc and t stability points even for broad states 
like e.g. the a and Kq mesons [33] . 



7. Masses of the 1 + four-quark states 



We reconsider the QSSR analysis in [22] using the previ- 
ous set of QCD parameters where the values of the gluon 
and four-quark condensates are about a factor 2 higher 
here. We shall introduce the renormalization group in- 
variant (i/'V') condensate [34] : 



Mi 



(logi//A) 



2/-/3i 



(25) 



with [311: A« = -(263± 7) MeV . We neglect the small 
Q^-dependence of asiip^) and {gipGip) as well as the 
anomalous dimension of the four-quark current, which 
should have small effect in the mass determination. The 
four-quark meson can be described by the diquark anti- 



6 



diquark operators 0: 

uICjf,dt{uaCd[ + UbCd^) + 



(26) 



where ... denotes an interchange between 7„ and 1 or 75. 
The QCD expression of the LSR of the corresponding 
correlator can be written to LO in as 1^: 



£4«(t) = 



where, for yyf^ 



dt [cof^ + C4t^ + CQt + cs) e 



-Cm + C12T 



(27) 



1 



Co 



C6 



ClO 



C12 = 



184327r6 

187r2 ' 

1927r2 
32 



C4 = 
5 



184327r6 



8647r2 



247r2 



81 S/DTT^ 



(28) 



from which one can deduce the ratio 7?.„(t) and Rn{tc) 
defined in Eqs. pH]) and pT|) used for extracting the 
meson mass. The result of the analysis is shown in Fig. 2] 
from which we can deduce in units of GeV: 



Miqi < 2.4(1) : LSR ; IVhqi = 1.70(4) 
for tc ~ (2.5 ~ 3) GeV^. 



FESR , (29) 



where one can notice that the LSR result does not have 
T- and ic-stabilitics such that only an upper bound can 
be extracted from the positivity of the spectral func- 
tion. This result is slightly higher than the one in [55] 
due to the consideration of the violation of factorization 
in the estimate of the high-dimension condensates here. 
We have repeated the analysis for the 77^^. We show 
in Fig [5^) the result of the analysis from the LSR mo- 
ment TZq{t) which presents r-stability for tc > 4: GeV'^ 
but no tc-stability. In the expanded form, FESR Ro{tc) 
presents a zero for 2.5 < tc < 7 GeV^ due to A^o(ic) 
signaling large non-perturbative effects for this quan- 
tity. This zero shows up like a minimum in tc in the 
non-expanded form of Ro(tc) given a mass M « 1.2 
GeV, but at too low value of ic ~ 2. GeV^, where the 
OPE does not converge H. At this ic-values, the LSR 
does not also show a r-stability. Therefore, we do not 
retain this result from our analysis. Ri{tc) moment has 
a much better behaviour but (unfortunately) does not 
show a tc-stability [see Fig. [Dd)]. The most conservative 
result from both LSR and FESR is: 



1.64 < M4g3 < 2.1(1) : > 4 GeV^ 



(30) 



*In principle, these operators should mix under renormalizations 
| 36| and one should built their renormalization group invariant 
(RGI) combination for describing the physical state. 
^In |22l . this minimum corresponds to a slightly higher value of 
the mass due to the different choice of the QCD set of parameters. 




Figure 5. Predictions of the four-quark state mass In GeV using 
the current In Eq. ( I3H : a) LSR: domain spanned by the upper bound 
versus r using the largest range spanned by the QCD parameters given 
after Eq. |[7j ( red: continuous curve): Mass using the central values 
of QCD parameters for tc — 3 GeV (green: dashed), — 4 GeV^ 
(black: dotted-dashed), ^ 5 GeV^ (blue: dashed), b) FESR versus 
tc using the central values of QCD parameters. 



where the upper bound comes from the positivity of 
the LSR moment. Currents with higher number of 7 
matrices have been also considered in |22j : 

^2m,4m = '>^lGYl5M<Gaf,^-f5d[±UbCcrfj,^'^5d^+...) , (31) 

where ... is an interchange between 7"^ and (7^,^. The 
analysis of the corresponding correlators lead to similar 
results (within the errors) than the ones in Eq. (|5S|) : 



M, 



M, 



and M, 



M, 



»)2f. 



(32) 



For definiteness, we consider that we have only one 4- 
quark state coupled to a RGI operator which results 
from a combination of these different operators expected 
to mix under renormalizations [36| . We fix the four- 
quark state mass to the value from Eq. ((29)l : 



= 1.70(4) GcV 



8. Masses of the 1 + molecules 



(33) 



The possibility to form two 1 + molecules with the two 
operators: 



(?A75-0)(V'757''V') , 



vvpa 



(S757prf)(d7"^) - (d757pw)(w7"rf) (34) 



has been discussed in [23], where the meson associated 
to the 1st (resp 2nd) operator is expected to decay into 
?77r, r(Tx (resp. p7r, 6i7r). Using the QCD expression 
given in [23] without the one-direct instanton contribu- 
tion which is largely affected by the uncertainty of the 
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3 i.i 4 4.3 5 Si. 5 6. .5.5 



Figure 6. Mass of the molecule state associated to the current Jt^ 
Eq. l|34j versus tc using the largest range spanned by the QCD 
parameters given after Eq. ([7j (red: continuous curve). 



instanton density p appearing as p^, we use the lowest 
FESR moment Ro{tc) for the analysis of the meson mass 
Mfnoii associated to the current where a ic-stability 
is obtained at — 5 GeV^ (Fig. [Q, while the LSR mo- 
ments 1^0,1,2 do not give conclusive results. The mass 
Mmoi2 associated to the current increases with tc 
for the FESR, while the LSR moments decrease slowly 
with the sum rule variable t but present a zero around 
r = 0.7 - 0.8 GeV"^ (not shown in Fig. [T]), while it 
presents ic-stability around 3 GeV^. The optimal re- 
sults are: 



moll 



1.41(3) GeV, 



3l2 



(1.2- 1.4)GeV,(35) 



which agree within the errors with the ones in . How- 
ever, we can consider that the two operators in Eq. ([34]) 
mix under renormalization such that only one physical 
state couples to the corresponding RGI operator, which 
we fix to be: 



M^oi = 1.3(1) GeV 



(36) 



9. Nature of the 7ri(1400), 7ri(1600) and 7ri(2015) 

From your previous analysis and for a further phenome- 
logical use, we shall consider one hybrid, one four-quark 
and one molecule states below 2 GeV, with the masses 
from Eqs. ([241), (m and (l36|) in GeV: 



Mh - 1.81(6) , Ahg = 1.70(4) , = 1.3(1) . (37) 

The 7ri(1400) is intriguing as it is seen to decay into rjn 
and rj'iT but not into pir, biir [1] . The later decays 
being also expected for an hybrid state ||9j. Looking 
at our bare (unmixed) mass predictions in Eq (j37p . one 
may expect that the molecule state is the most probable 
candidate. We consider that the suppression of the pw 
decay can be due to a mixing of this molecular state with 




0.3 0.5 1 1.1 1.2 l.i 1.4 



Figure 7. Mass of the molecule state associated to the current 
in Eq. ( I34t from LSR versus the sum rule variable t= 0.1 to 0.6 
GeV^^ and 0.8 to 1.4 GeV~^ and for different values of tc using the 
central values of the QCD parameters given after Eq. |7J and for two 
LSR moments TZq (continuous curve) and TZi (dotted-dashed curve): 
tc = 2.5 GcV^ (green), = 3 GoV^ ( red) and = 3.5 GcV^ (blue). 



the four-quark or/and hybrid states. We use a minimal 
two-component mixing: 



TTi (1400) = cos 9mol\mol) + sin 9rnol\X) 

7ri(1600) = - sin dmoi\mol) + COS 9moi\X) 



(38) 



where X is a four-quark or hybrid state. The "best fit" 
is obtained for the sets : 

{MmouMx) - (1.2- 1.3 ; 1.70,1.74) GeV ^ 

Braol ^ -(11.7 ± 2.2)° , (39) 

which is slightly favours a molecule/four-quark mixing. 
This result may suggest that the 7ri(2015), quoted in 
the extended version of PDG [T], can have more hybrid 
state in its wave function. One also expects that the 
non-exotic meson 1 as well as the isospin partners of 
these 1 exotics should be almost degenerate in masses 
with these exotic states, while some of their radial exci- 
tations have masses of the order of optimal continuum 
threshold: 



M' w V^I^i (1.7- 2.2) GeV 



(40) 



Further experimental ad theoretical tests of this mixing 
scheme are required. 



°For a review on different experimental results and related prob- 
lems, see e.g. 1371. 
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